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ABSTRACT 

We report the detection of a large mass planet orbiting around the KO metal- 
rich subgiant HD38801 {V — 8.26) by precise radial velocity (RV) measurements 
from the Subaru Telescope and the Keck Telescope. The star has a mass of 
I.36M0 and metallicity of [Fe/H]= +0.26. The RV variations are consistent 
with a circular orbit with a period of 696.0 days and a velocity semiampli- 
tude of 200.0ms~^, which yield a minimum-mass for the companion of lO.TMjup 
and semimajor axis of 1.71 AU. Such super-massive objects with very low- 
eccentricities and hundreds of days period are uncommon among the ensemble 
of known exoplanets. 

Subject headings: stars: individual:HD38801, stars: planetary system - techniques: 
radial velocities - techniques: spectroscopy 
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Introduction 



Since 1995, precise RV measurements have unveiled m ore than 400 extrasola r planets 
from surveys that include roughly 3000 nearby stars (e.g., 



Udry fc Santos 



20071 ). The 



planets show a great diversity in their physical attributes: masses of the planetary 
companions range from 4M^{r^ O.OlSMjup) to larger than 13Mjup, with semimajor axes 
of 0.02AU - 6AU and orbital eccentricities of - 0.9. A large number of planets have 
enabled us to discuss correlations between each orbital element and their masses in terms 
of planet formation and orbital evolution. For example, the diagram of orbital eccentricity 
against semimajor axis shows that eccentricities of planets with semimajor axis larger than 
0.1 AU are almost uniformly distributed from to 1, suggesting orbital evolution due to 
gravitational interaction between planets, while those inside 0.1 AU tend to be damped 
below 0.1 due to tidal circularization. The diagram of eccentricity against planet mass also 
shows another possible hint for planet formation that super- massive planets (> 5Mjup) 
tend to have relatively high eccentricity (> 0.3), while less massive planets have a wide 
range of eccentricities. 

As the number of planet surveys has enlarged, diversity of planet-host stars has also 



been explore d. Several surveys have focused on specified ste 



stellar mass (iBonfils et al. 



and metallicity ( iFischer et al 



2005 



Apps et al. 



2005 



2010 



Sozzetti et al. 



Sato et al. 



ar properties, particularly 



2008 



Johnson et al. 



2003) 



20091 ) ■ in order to investigate correlations 



between the stellar properties and planetary parameters such as stellar mass and occurrence 
rate of planets. As a result from these surveys, the statistical correlations between stellar 
properties and planet parameters have begun to be unveiled. For example, occurrence rate 
of Jovian planets increases with st ellar metallicity fro m 3% for stars with [Fe/H] < to 



25% for those with [Fe/H] > +0.3 flFischer et al. 



20051 ) . and also increases with stellar mass 



from 2% around M dwarfs (< 0.6Mq) to approximately 9% around F and A stars (1.2 - 1.9 
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metallicity and mass ( iLovis fc Mayor 



M©; I Johnson et al.ll2007|). Planetary mass could also be positively correlated with stellar 



20071 ). Impacts of stellar properties on characteristics 



of planets, such as orbital eccentricity and multiplicity, are also expected to emerge as the 
number of planets grows. 

The N2K program is a large scale international exoplanet-search project started in 
2004, which consists of a sample of 2000 metal-rich solar-type stars (IFischer &: Valenti 
20051 ). The search originally targeted short-period planets with a high-cadence observational 
strategy. However, it has also detected intermediate-period (18 — 3810 days) ones thanks 
to the long-term observations. From the collective N2K surveys, we have discovered 
22 planets so far, which have a wide variety of mass (0.22 - 13.1 Mjup) and orbital 
parameters (a = 0.04 - 4.9AU, e = - 0.7), a nd orbit around a variety of host stars 



(-0.11 < 



Peek et al 



e/Hl < 0.37,0.88 < M/Mq < 1.31; 



Wright et al 



2007 



Johnson et al. 



2006 



20091 ). The planets discovered from our surveys can thus contribute to our 



general understanding of planet formation and evolution depending on stellar properties. 

Here we report the detection of an exoplanet orbiting the metal-rich solar type subgiant 
HD38801 discovered from the N2K sample at the Subaru Telescope and Keck observatory. 
The planet has a minimum mass of lOMjup and a low eccentricity. Such a planet has 
been rarely discovered so far around solar-type dwarfs and subgiants. We describe the 
characteristics of HD38801 in ||2l In §3l we present our observations, and the orbit of 
HD38801b. We summarize our results in ^ and discuss formation scenarios of such a 
low-eccentricity super- massive planet. 
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Stellar characteristics 



HD38801 (HIP 27384) is listed in the Hipparcos catalogue (ESA 1997) as a KOIV star 
with a visual magnitude of V = 8.26, and a color index B — V = 0.873. The Hipparcos 
parallax of 10.06 mas corresponds to a distance of 99.4 pc which yields the absolute 
visual magnitude of My = 3.27. A high- resolution spectroscopic analysis described in 
Valenti fc Fischerl (120051 ) derived an effective temperature Te// = 5222 ± 44 K, a surface 



gravity logg = 3.84 ± 0.1, rotational velocity vsini = 0.54 ± 0.5km s ^, and metallicity 
[Fe/H] = 0.26 ± 0.03 dex for the star. The bolometric luminosity of L,^tnr = 4. 56 ± 0.47^0 



19961 ). From the 



was derived from the My and a bolometric correction of —0.195 (iFlowerl 
bolometric luminosity, parallax and T^jf, we derived a stellar radius of 2.53 ± 0.13Rq based 
on the Stefan-Boltzmann relation. 

In order to estimate a stellar mass, we interpolated the m etallicity, effective 



temperature, and luminosity onto the stellar interior model grids flGirardi et al 



Johnson et al. 



2002h. We 



fl2007h 



adopted the three-dimensional interpolation method described by 
and applied 7% uncertainty based on a comparison among different stellar model grids. 
Using the SME-derived stellar parameters, we estimated mass of 1.36 ± O.O9M0 and an age 
of 4.67 ± 2.56 Gyr for IID38801. We measured Shk of —0.174, which is core emission in 
the Ca II HK lines relative to the continuum, to obtain chromospheric activity. We derived 
logR'fjj^ = —5.0, the ratio of flux from Shk to the bolometric stellar flux, which indicates 
that the star is chromospherically inactive. However, it is known that subgiants show 
"jitter" of 4-6 ms~^, which is velocity scatter in excess of internal errors due to astrophysical 
sources such as pulsat ion and rotational ra odulation of surface features, even with such a 



low value of log R'hk ( j Johnson et al. 



r2 

error 



which is added in quadrature (a, 
HD38801 when we fit a Keplerian orbit to the RV data. 



20071 ). We therefore adopt a jitter value of 6 ms ^ 

jitter J to the formal velocity uncertainties for 
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3. Doppler observations and Keplerian-fit 



The Subaru Telescope's component of the N2K survey uses the High Dispersion 



2OO2I). We adopted 



Spectrograph (HDS) on the 8.2-m Subaru Telescope ( iNoguchi et al. 
the setup of Stdl2b, which simultaneously covers a wavelength region of 35OO-6IOOA by 
a mosaic of two CCDs. The slit width was set to 0.8" for the first 3 observations (2005) 
and 0.6" for the last 15 data, giving a reciprocal resolution (A/AA) of 45000 and 55000, 
respectively. Using this setup, we can obtain a signal-to-noise ratio of S/N ~ 150 pixeP^ at 
5500 A with 110-300 sec exposure for our typical targets, V ~ 8.5 stars. We used an iodine 



absorption ce 



( iKambe et al 



to provide a fiducial wavelength reference for precise RV measurements 



2002 



Sato et al 



2OO2I ). Although the cell was originally installed just behind 



the ent rance of the spect rograph, it was moved in front of it in 2006. We used an analysis 



code of 



Sato et al. 



( I2OO2I ) for modeling iodine superp osed spectra and de rivin g stellar RVs, 



which is based on the standard method developed by 



Butler et al. 



(11996h and 



Valenti et al. 



( 119951 ). Using the code and above spectrograph setting, we can achieve a Doppler precision 



of 4-5 ms ^. 



Between Sep 2006 and Jan 2009, we obtained 10 RV observations from the Keck 10-m 



telescope with the HIRES spectrograph (IVogt et al. 



19941 ). An iodine cell was placed in the 



optical path to provide both the precis e wavelength solution and the PSF 



Function) of spectrograph instruments ( iMarcy fc Butler 



1992 



Butler et al 



Poin t Spread 



19961 ). We used 



the B5 decker corresponding to a slit width of 0.86" and a wavelength resolution of 65,000. 
Appropriate exposure times between 150 and 240 seconds were taken and this results the 
S/N of 150 for respective observations. The typical precision of recent Keck observation is 
1.5 ms"'^. 

We made the first observations for HD38801 at the Subaru Telescope and detected 
large RV variations. After that, the star has been monitored at the Subaru Telescope and 
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the Keck Telescope for 3 years (from November 2005 to November 2008) and a total of 20 
data points were obtained. The observation dates, RVs, and instrumental uncertainties are 
listed in Table [2] and plotted in Figure [TJ 



To derive the best -fit Keplerian or 



algorithm (AMOEBA; 



jit to the RV data, we adopted the down-hill simplex 



Nelder fc Mead 



19651 ) to reduce which has multiple parameters 
(orbital period P, time of periastron passage Tc, eccentricity e, velocity semiamplitude Ki, 
and argument of periastron u) and obtained the best fit parameters of single Keplerian orbit. 
The stellar jitter of 6 ms^^ was added in quadrature to the measurement uncertainties. At 
first fit, when all of these parameters were let free, we obtained very low eccentricity of 
0.04 and poorly constrained uj. So we made the another fitting with the fixed values of e 
and UJ to 0. The best fit parameters thus obtained are listed in Table [3l and the Keplerian 
fit is overplotted on the RV data in Figure [H We obtained a period of P = 696.3 ± 2.7 
d and a velocity semiamplitude of 200.0 ± 3.9 ms~^. The rms scatter to this fit is 6.5 
ms~^ and (x^)^^^ = 1-2 as a reduced x- The uncertainties in the orbital parameters were 
estimated by a bootstrap Monte Carlo approach. The RV residuals from best-fit Keplerian 
curve are scrambled and added back to the original measurements, and then re-fit and 
derive new parameters. Iterating 10,000 times, we obtained la value of each parameters. 
With a stellar mass of 1.36 ± O.O9M0, we obtained a minimum-mass for the companion of 
Msinz = 10.7 ± 0.5Mjup and semimajor axis of = 1.70 ± 0.03AU. 

As seen in Figure [H the residuals to the Keplerian fit in JD24543 00-2454500 ex hibit 



19821 ) to the 



larger scatter than other parts. We performed a periodogram analysis (iScargld 
residuals. However, we found no significant peaks in the periodogram at this stage. More 
dense sampling of data will help verify the existence of any periodicity in the residuals and 
its origin. 
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4. Summary and discussion 

We here reported the detection of a large mass planet in an almost circular orbit 
around a high metallicity KO IV type star HD38801 from the precise RV observations at the 
Subaru Telescope and the Keck Telescope. The star has a mass of 1.36Mq and metallicity 
of [Fe/H]=0.26. The planet is a new sample around a metal-rich and relatively high-mass 
star, which makes an important role to investigate dependence of planets on host star's 
properties. 

One remarkable feature of this planet is its low orbital eccentricity (e ~ 0) despite of 
the large mass (Msini = 10.7Mjup) and the intermediate orbital distance (a = 1.7AU). 
Such a planet with low eccentricity (e < 0.1) and large-mass (Msini > lOMjup) has 
been rarely discovered so far in the regions of 0.1 AU < a < 5 AU around FGK dwarfs 
and subgiants. Only HD16760b (Msini = 13Mjup,a = 1.08 AU, e = 0.084), which w as 



recently discover ed from the Subaru Telescope and the Keck Telescope (IBouchy et al. 



2009 



Sato et al 



20091 ). falls into this category except for HD38801b. In the core-accretion model, 
massive companions are generally formed with circular orbits because of protoplanetary 
disk's circular motion. The circular orbit of IID38801b is consistent with this scenario. 
Tidal interaction between the host star and the planet cannot explain the circular orbit 
because the stellar radius is only 2.5Rq (~ O.OIAU), which is too small to circularize the 
planetary orbit at 1.7 AU. Contrary to HD 38801b, most of the intermediate-period (> 18 
days) and large-mass (> 5Mjup) planets discovered around FGK dwarfs and subgiants 
so far tend to reside in eccentric orbits. The origin of such planets still remains to be 
solved. It is generally thought that eccentric planets are formed by planet-planet scattering 
in multiple planetary systems. The scenario suggests formation of multi super-massive 
planets in a system and predicts existence of scattered outer companions as massive as 
inner ones, which can be good targets for direct imaging. Other scenarios such as giant 
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impact ( iBaraffe et al.l 120081 ), and disk in stabilitv (lBossl Il998l ) and Kozai-mechanism when 



they are in binary systems (IKozai 



19621) have also been proposed to explain for formation of 



such massive eccentric planets. Interestingly, planets with similar properties to HD38801b 
have been discovered around intermediate-mass giants. Comparing populations of such 
planets between different types of stars would provide a hint on formation scenario and 
evolution for the planets. 

The high-metallicity and the evolutionary stage of subgiant for HD38801 is another 
interesting feature of this system. It is well known that the dete ction rate of giant plane ts 



20051) 



shows positive correlation with the metallicity of the host stars ( iFischer fc Valenti 
Two scenarios have been proposed to explain the origin of the correlation; one is that 
more metal-rich stars tend to form more gia nt planets and t he other is that the stellar 



surface is polluted by planetary debris (e.g., 



Ida fc Lin 



2004 - At the early stage of the 



main sequence, a stellar surface could be polluted due to the accretion of planetesimals 
and protoplanets from debris disk and metallicty of the stellar surface could be enhanced. 
However, if the star has deep convective envelope, the metal-enhanced surface can be mixed 
with original stellar material and diluted by the convective flow. Thus, the existence of a 
giant planet around metal- rich evolved star as the subgiant IID38801, which probably has 
deep convective envelop e, favors the former scen ario, that is giant planets tend to form in 



metal-rich environment (IFischer &: Valenti 



20051 ). It should be noted, however, that many 



planets have been discovered even in metal-poor ([Fe/H] < 0) giants, more evolved stars 
than subgiants , and no significant planet-metallicity correlation can be seen among them. 



Pasquini et al. 



(120071 ) proposed that the lack of metal-rich tendency is due to dilution by 
deeper convective envelope of giants than that of dwarfs and subgiants. To make sure the 
origin of the high metallicity of planet-harboring dwarfs, it is important to increase the 
number of planets discovered around evolved stars with < I.5M0, for which we can detect 
planets in all of the three evolutionary stages, dwarfs, subgiants, and giants, by precise 
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Doppler technique. 
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Table 1. Stellar Parameters for HD 38801 



Parameter Value(error) 



V 


8.26 


Mv 


3.27 


B-V 


0.873 


Spectral type 


KO IV(subgiant) 


Distance(pc) 


99.4 


re//(K) 


5222(44) 


logs 


3.84(0.1) 


[Fe/H] 


0.26(0.03) 


V sini(kms^^) 


0.54(0.5) 


M*(M0) 


1.36(0.09) 




2.53(0.13) 




4.56(0.47) 


age(Gyr) 


4.67(2.56) 


Shk 


-0.174 




-5.01 
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Table 2. RVs for HD38801. 



Observation Date RV Uncertainties 

(JD -2, 450, 000) (ms-i) (ms-i) Observatory 



3695.0902 


-159.0 


4.2 


Subaru 


3695.9883 


-151.3 


4.0 


Subaru 


3727.0304 


-115.4 


4.8 


Subaru 


3982.1185 


201.4 


0.9 


Keck 


3983.1122 


196.8 


0.9 


Keck 


3984.1088 


202.8 


0.9 


Keck 


4081.9436 


105.9 


4.0 


Subaru 


4082.9788 


101.9 


4.7 


Subaru 


4085.9284 


92.2 


4.8 


Subaru 


4087.0406 


89.0 


4.1 


Subaru 


4130.8895 


13.4 


1.2 


Keck 


4138.9172 


-6.9 


1.2 


Keck 


4344.1185 


-180.3 


0.8 


Keck 


4398.0809 


-154.4 


0.9 


Keck 


4427.9607 


-93.3 


1.1 


Keck 


4470.9194 


-19.9 


4.6 


Subaru 


4492.7365 


-3.1 


1.0 


Keck 


4523.7729 


59.0 


3.5 


Subaru 


4744.0407 


145.0 


4.1 


Subaru 


4787.0288 


89.6 


3.9 


Subaru 


4838.9756 


4.0 


1.2 


Keck 
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Table 3. Spectroscopic Orbital Solution for HD38801b 



Pcirameter Value(error) 



P(d) 

Tc(JD) 

Eccentricity 

^^(ms-l) 

ai sini(10-3AU) 

ap(AV) 

Msini(Mjup) 

Nobs 

rms(ms~^) 
Reduced(x^)^''^ 



696.3(2.7) 
2,453, 966.0(2.1) 
O(fixed) 
O(fixed) 
200.0(3.9) 
13(0.27) 
1.70(0.037) 
10.7(0.5) 
21 
6.5 
1.2 
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Fig. 1. — The RVs for HD38801 derived from the Subaru Telescope (/i/Zeci circles) and the 
Keck Telescope( openec? triangles). With an orbital period of 696 d, velocity semiamplitude 
of 200ms~^, and stellar mass of 1.36Mq, the indicating planet mass is Msini = 10.7Mj, 
the semimajor axis is 1.7AU. The eccentricity is fixed to 0. The rms to the fit is 6.5ms~^. 
The offsets of — 58.8ms~^ and 20.5ms~^ were applied to the Subaru RVs and the Keck RVs 
respectively. 
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Fig. 2. — The correlation diagram of the eccentricity and the minimum mass of the planets 
with above 14 days orbit. HD38801b is plotted with filled triangle. HD16760b(/j//e(i circle), 
main sequence stars(o]5enec? circles), and subgiants(o])ene(i triangles), and giantsf croggea) are 



also plotted in. The companion's properties of HD 16760 is referred from 



Sato et al. 



(|2009|) 
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